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PREFACE 


As  one  element  of  a  continuing  program  to  study  AGILE  communication 
requirements,  Defense  Research  Corporation  has  been  analyzing  the  trans¬ 
mission  of  radio  signals  through  tieavily  vegetated  (jungle)  terrains. 

The  major  objective  of  this  work  is  to  devise  theoretical  models  that  can 
be  used  to  predict  the  signal  strengths  to  be  expected  in  a  variety  of 
practical  situations.  Our  prograts  to  accomplish  this  examines  the 
electromagnetic  characteristics  of  particular  models  of  the  propagation 
medium  and  then  tests  the  usefulness  of  the  model  by  its  ability  to 
represent  experimentally  available  results  in  a  meaningful  way. 

Beginning  with  the  simplest  possible  model,  we  are  introducing 
increasing  levels  of  sophistication  as  required  by  the  comparison  with 

experiment.  For  example,  in  the  initial  phase,  the  transmission  medium 

« 

was  represented  as  an  infinite,  homogeneous  region,  (see  Ref.  2  of  the 
present  report),  characterized  only  by  a  complex  dielectric  constant.  In 
this  report,  a  somewhat  more  complicated  model  is  analyzed:  a  uniform 
slab,  of  thickness  corresponding  to  the  height  of  the  jungle,  bounded  on 
one  side  by  an  Infinite  lossless  region  (air),  and  on  the  other  side  by 
an  infinite,  lossy  region  (ground).  Further  refinements  will  be  intro¬ 
duced  in  the  future. 

This  program  is  under  the  direction  of  Dr,  B.  A.  Lippmann.  Tlie  work 
reported  on  here  was  carried  out  by  Drs.  D.  L.  Sachs  and  P.  J.  Wyatt,  with 
substantial  assistance  from  Dr.  P.  J.  Redmond, 
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ABSTRACT 


A  theoretical  determination  of  the  path  loss  In  radio  propagation 
through  jungle  Is  obtained  by  considering  the  jungle  as  a  homogeneous 
conducting  dielectric  slab  on  a  flat  earth.  An  exact  Integral  is  ob¬ 
tained  for  the  vertical  component  of  the  electric  field  within  the 
jungle  due  to  a  vertical  electric  dipole  within  the  jungle.  An  analytic 
evaluation  of  the  integral  leads  to  an  approximate  expression  for  the 
field  and  aa  estimate  of  the  error.  When  this  accuracy  is  insufficient, 
a  numerical  evaluation  of  the  integral  is  performed.  The  combination 
of  analytic  and  numerical  techniques  leads  to  the  evaluation  of  the 
electric  field  and  thereby  path  loss  to  any  accuracy  desired.  Experi¬ 
mental  measurements  of  path  loss  in  a  Thailand  jungle  are  compared  with 
the  calculations.  The  results  agree  within  a  standard  deviation  of  6  db 
when  the  jungle  conductivity  is  taken  as  0.15  mllllmho/m. 
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I.  INTRODUCTION 

The  practical  difficulties  associated  with  radio  conaunlcatlons 

in  regions  of  dense  vegetation  have  been  recognized  for  some  tlme.^ 

The  vegetation  has  a  high  water  content  so  that  the  jungle  medl<mi  Is 

lossy.  If  the  radio  waves  propagated  along  a  direct  path  through  the 

2 

oMdium  thjB  ricelved  signal  would  be  exponentially  attenuated.  With 
reasonable  values  of  the  conductivity  (a  =  100  ymho/m)  the  attenuation 
is  approximately  0.2  db/m  at  one  megacycle  and  is  even  greater  at 
higher  frequencies.  Such  a  high  attenuation  would  effectively  prevent 
radio  conmunlcatlon  over  distances  greater  than  a  few  thousand  feet. 
Since  radio  communication  docs  exist  over  greater  distances  than  this» 
albeit  quite  poorly  on  occasion,  there  has  been  speculation  that  the 
radiation  propagates  for  the  most  part  in  the  air  above  the  jungle. 

This  paper  derives  from  first  principles  expressions  for  the  pro¬ 
pagation  between  two  antennas  immersed  in  a  jungle-like  medium.  The 

theory  is  compared  with  experiments  conducted  in  a  Thai  jungle  by 
3  4 

Jansky  and  Bailey  *  and  the  apparent  anomaly  referred  to  appears  to  be 
explained. 

The  electrical  properties  of  a  jungle  are  extremely  complicated. 
Most  of  the  space  is  occupied  by  air  with  properties  close  to  the 
vacuum.  The  trees  and  other  vegetation  have  a  higher  conductivity  and 
dielectric  constant  and  are  distributed  more  or  less  randomly  over  the 
surface.  Fortunately,  at  long  wavelengths  the  propagation  character¬ 
istics  of  an  electromagnetic  wave  are  not  sensitive  to  this  fine  detail. 
An  electroouignetlc  wave  propagating  through  a  medium  is  only  sensitive 
to  some  average  properties  of  the  medltm  where  the  average  Is  taken 

2 

over  regions  whose  linear  dimensions  are  of  the  order  of  a  wavelength. 

We  therefore,  consider  an  approximation  to  the  jungle  consisting 
'  of  a  tinlform  slab  of  fixed  height  with  permitlvlty  and  conductivity 
d.  bounded  by  a  flat  earth  surface  with  permitlvlty  e  and  conductivity 

J  €9  - 

3  and  above  by  the  air  with  vacuum  properties.  This  picture  is 
g 

reasonable  provided: 


m 


1.  The  flifctuatlon  la  the  number  of  trees,  etc..  In  an  area 
one  waveli^ngth  squared  is  small  compared  to  the  total  number 
of  trees  in  this  area. 

2.  If  the  height  of  the  - ungle  is'  larger  than  a  wavelength  it 
is  necessary  that,  within  the  jungle,  the  average  electrical 
properties  do  not  vary  significantly  with  height. 

3.  The  transition  region  between  the  air  above  and  the  jungle 
must  be  small  compared  to  a  wavelength. 

These  conditions  are  prob/bly  sufficient  to  guarantee  the  quali¬ 
tative  validity  of  the  theory.  For  quantitative  validity  it  would  also 
be  necessary  that  the  terrain  be  reasonably  level  over  the  range  and 
that  the  properties  of  the  jungle  be  reasonably  uniform  over  the  range; 

In  this  study  it  was  Intended  to  determine  the  range  of  validity 
of  this  idealized  slab  model  by  comparing  the  results  of  calculations 
with  experlmenf-al  data.  For  the  wavelength  region  studied  so  far  (wave¬ 
length  from  3  to  50  meters)  we  have  obtained  satisfactory  agreement 
between  theory  and  experiment. 

The  problem  of  the  propagation  of  electromagnetic  energy  through 
a  many- layered  medium  has  been  studied  extensively. For  this 
general  class  of  problems  it  is  relatively  easy  to  obtain  a  formal 
solution  to  the  problem  in  which  the  field  strength  at  the  receiver  is 
represented  by  an  Integral.  We  have  obtained  such  an  integral  represen¬ 
tation  for  the  particular  geometry  of  our  problem. 

Unfortunately,  the  Integral  representation  of  the  solution  is 
sufficiently  complicated  that  considerable  effort  is  required  before 
numerical  results  can  be  obtained  in  a  form  suitable  for  comparison  with 
experiment.  We  have  performed  these  calculations  with  sufficient  pre¬ 
cision  so  that  any  discrepancy  between  the  mathematical  results  and  the 
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experimental  data  is  due  either  to  inadequacies  in  the  physical  model 
or  uncertainties  In  the  experiment. 

At  larger  ranges  it  is  possible  to  obtain  an  approximate  analytical 
evaluation  of  the  Integral  and  an  estimate  of  the  error  involved  in  the 
approximation.  This  acyozptotic  solution  corresponds  to  a  lateral  wave 
propagation  mode,  or  "treetop"  mode.  This  mode  may  be  described  in 
terms  of  a  wave  packet  (ray)  which  propagates  up  to  the  treetops,  strik¬ 
ing  the  boundary  near  the  critical  angle  for  total  internal  reflection. 

It  then  propagates  in  the  air  until  it  reenters  at  an  angle  near  the 
critical  angle.  The  field  is  attenuated  exponentially  only  along  a 
path  from  the  transmitter  to  the  treetops  and  from  the  treetops  to  the 

2 

receiver.  While  propagating  through  the  air  the  field  falls  off  as  1/r  . 

This  asymptotic  formula  la  capable  of  explaining  most  of  the  experi¬ 
mental  data.  At  the  smaller  ranges  the  error  becomes  too  large.  For 
this  region  the  original  Integral  was  evaluated  by  numerical  quadrature. 

.\8  a  check  on  the  consistency  of  the  two  procedures  the  numerical 
quadrature  was  carried  up  to  ranges  where  the  asymptotic  formula  was 
valid.  In  the  region  of  overlap  the  two  methods  gave  identical  results. 
Examples  where  all  of  these  numerical  techniques  are  required  are  given 
in  Sec.  Ill,  and  the  transition  regions  are  discussed. 

In  Sec.  II  the  formal  mathematical  solution  is  presented  and 
it  is  demonstrated  that  this  solution  is  consistent  with  our  description 
of  the  lateral  waves. 

In  Sec.  Ill  the  path  loss  for  a  variety  of  experimental  data  is 
compared  with  the  theoretical  predictions.  To  date  we  have  considered 
only  the  cases  where  both  the  receiving  and  transmitting  antennas  were 
in  the  jungle  medium.  There  is  a  considerable  amount  of  data  for  this 
case  and  it  is  in  this  situation  where  the  greatest  difficulties  in 
communication  might  be  anticipated.  The  calculations  considered  only 
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the  case  of  vertical  polarization,  a  frequency  range  from  6  to  100  Mc/sec, 
and  a  range  up  to  one  mile.  Within  these  restrictions  the  available 
data  provided  by  Jansky  and  Bailey  was  considered.  The  theoretical 
results  depend  most  critically  on  the  jungle  conductivity  and  this  was 
the  only  parameter  varied.  The  theoretical  and  experimental  results 
agree  to  within  a  standard  deviation  of  6  db  when  the  jungle  conduc¬ 
tivity  is  taken  as  0.15  milllmho/m.  This  agreement  is  considered  satis¬ 
factory.  Jansky  and  Bailey  performed  Independent  measurements  corres¬ 
ponding  to  the  same  range,  height  of  antennas,  and  frequency  with  results 
varying  by  as  much  as  13  db. 

The  average  conductivity  of  the  Thailand  jungle  has  not  been 
measured,  and  there  are  very  few  measurements  of  this  kind  available 
for  any  jungle.  The  value  we  have  obtained  is  reasonable  and  is  con- 

g 

slstent  with  experimental  results  obtained  in  California. 

The  slab  model  predictions  appear  to  provide  a  satisfactory 
correlation  with  the  experimental  data.  Indeed  the  model  works  well 
for  wavelengths  shorter  than  we  had  at  first  anticipated  and  we  intend 
to  extend  our  calculations  to  still  shorter  wavelengths  in  order  to 
determine  when  and  if  the  model  seriously  breaks  down.  We  also  Intend 
to  consider  horizontal  polarizations  and  cases  where  one  or  both  of  the 
antennas  are  above  the  jungle. 

The  details  of  the  numerical  evaluations  and  the  derivation  of 
the  formal  solution  are  contained  in  appendices. 

II.  PHYSICS  OF  THE  SLAB  MODEL 

The  mathematical  problem  posed  by  the  idealized  physical  situation 
possesses  an  exact  formal  solution.  It  suffices  to  consider  only  the 
formula  for  the  z-component  of  the  electric  field,  since  the  other 
field  components  are  readily  obtained  from  E^.  The  z  component  of  the 
electric  field  at  the  receiver  is  given  by 
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E  ■  P  /  H  (oR)  F(a)  da  volt/m 

^  ^a,  ° 

where  P  ■  (power  radiated  in  kw)^^^  and  f  is  the  frequency  in  Me. 

n 

For  both  the  receiving  and  transmitting  antenna  in  the  jungle 

.  r  -X  X  /  -X  Y  -X  (2H-Y)  -X  (2H+X) 

F(o)  -  ~  e  •^+|veJ+VeJ  H-VVe^ 

Xj  I  \  g  a  a  g 

1/ 


(1) 


+  V  V  e 
a  g 


-Xj  (2H-X)\ 


(2) 


where 


V 

Xj  +  nj  x^ 


2  2 
n  X.  -  n.  X 

g  1 _ L_g. 

2^2 
n  X.  +  n.  X 

g  j  j  g 


where  n  *  e  +  18  lo  (mllliffiho/m)/£(Mc).  (See  Fig.  1.)  This  Integral 
represents  the  solution  in  the  jungle  as  a  superposition  of  cylindrical 
waves  which  are  successively  reflected  from  the  jungle-air  and  jungle- 
ground  interfaces  with  reflection  coefficients  and  respectively. 
The  variable  of  integration,  a,  is  the  cylindrical  propagation  constant. 


Distances  are  measured  in  units  c/w  «  150/ (irf)  »  1/K 

o 

and  R«Kr,  H-Kh,  X"K  |z-z  | ,  Y  »  K  (z+z  ),  Z-K  z,Z  «K  z  . 

o*  o*  o'  o'^  o'  o^^  o’ooo 

The  attenuation  factors  for  different  paths  are  described  in  terms  of 
the  variables 


and  the  square  roots  are  defined  so  that  the  real  parts  of  the  x's  are 
positive. 
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The  various  terms  In  Eq.  (2)  have  a  simple  Interpretation.  The 

3C  • 

first  term  e  J  represents  the  field  produced  by  the  dipole  In  the 

-X4Y 

absence  of  the  Interfaces.  The  next  term  proportional  to  V^e  ■' 
represents  those  plane  waves  which  are  reflected  an  odd  number  of  times 
with  the  first  reflection  occurring  at  the  ground.  The  remaining  terms 
represent  cases  where  there  Is  an  odd  number  of  reflections  with  the 
first  reflection  at  the  air  Interface,  an  even  number  of  reflections 
with  the  first  reflection  at  the  air  Interface,  and  an  even  number  of 
reflections  with  the  first  reflection  at  the  ground  Interface.  The 
numerators  represent  the  results  of  the  first  one  or  two  reflections 
respectively  and  the  denominator  sums  the  contributions  from  all  sub¬ 
sequent  pairs  of  reflections  from  the  ground  and  air  Interfaces. 


For  most  values  of  the  parameters  the  Integrand  In  Eq.  (1)  Is 
strongly  peaked  at  a  =  1  corresponding  to  a  ray,  called  the  lateral 
wave,  propagating  at  the  angle  of  total  Internal  reflection.  In  such 
cases  the  Integral  may  be  approximated  by  a  simpler  Integral  describing 
this  ray.  The  new  Integral  can  be  evaluated  easily  by  a  Gaussian 
quadrature.  approximate  analytical  evaluation  of  the  new  Integral 
can  also  be  performed.  Such  solutions  are  presented  la  Appendix  III 
and  the  domains  of  their  validity  are  discussed. 


Uhen  the  methods  described  above  break  down,  It  becomes  necessary 
to  evaluate  the  original  Integral,  Eq.  (1),  numerically.  This  Is  a 
difficult  task  because  the  Integrand  Is  oscillatory.  Methods  for 
accurately  evaluating  such  oscillatory  Integrals  are  available  and 
they  were  applied  to  this  type  of  problem  for  the  first  time  In  this 
study. 


When  the  solution  can  be  expressed  In  terms  of  the  superposition 
of  a  set  of  plane  waves  with  propagation  vectors  lying  In  a  small  region 
of  k-space  It  Is  most  appropriate  to  Interpret  the  solution  In  terms  of 
a  wave  packet  which  suffers  successive  reflections.  In  order  to  gain 
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some  Insight  Into  the  nature  of  these  solutions  we  consider  the  behavior 
of  a  wave  packet  after  a  single  reflection. 


Consider  a  wave  packet  with  average  propagation  vector  k  Incident 
on  a  plane  boundary  as  in  Fig.  2.  For  simplicity  we  assume  the  fields 
are  independent  of  y.  The  incident  wave  packet  can  then  be  described  by 

E  ^’^^(x,z,t)  =  /  dk  dk  A(k  ,k  )  expi(k  x  +  k  z  -  tot)  (3) 

z'**'  xz  x*z  *^  x  2 


where  to  is  a  function  of  k^  and  k^  and  the  function  A  is  strongly  peaked 
at  ^  If  we  expand  the  phase  factor  about  the  point  k  *  we 

obtain 


„  me ,  ^ . 

(x,z,t) 


ik°x+  Ik^^z  -  ito°t 
X  z 


(4) 


where  3(o/3k  and  3<o/3k  are  components  of  the  group  velocity  evaluated 
at  k°  and 


APPARENT  PLANE  OF  REFLECTION 


Figure  2.  Displacement  of  Wave  Packet 
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The  reflected  wave  Is  given  by 


E*^®^(x,z,t)  •Jdk^dk^  A(k^,k^)  lv(k^,k^)  |exp  i  (k^x-k^z-(Dt+i|i(k^,k^))  (6) 

where  the  reflection  coefficient  is  jvle^*^.  If  we  assume  that  the 
amplitude  of  the  reflection  coefficient  varies  sufficiently  slowly  so 
that  it  may  be  taken  outside  the  Integral  we  obtain 


,  l(k  °x  -  k  °z  -  a)°t  +  (j)(k  ®,k  °)) 

„ref,  o.Ovi  X  z  ^  x’z  ' 

E  (x.z.t)  -  |V(k  ,k  )|  e 

X  » 


.  f/v  -  —  t  -f  -  z  -  t  + 

3k  ^  3k  *  3k  '  ^  3k  / 

'XX  Z  Z' 


where  f  is  the  same  function  that  occurs  in  Eq.  (4)  and  the  derivatives 
are  again  evaluated  at  . 

Consider  now  an  incident  wave  packet  which  first  strikes  the  boun¬ 
dary  at  t  «  0  and  at  this  time  is  nonvanishing  only  for  negative  x  and  z. 
These  conditions  require 


f(a,3)  ■  0  if  a  >  0  or  B  >  0. 

By  examining  Eq.  (7)  we  see  that  the  reflected  signal  will  not  appear 
in  region  I  (Illustrated  in  Fig.  2),  until  a  s\ifflclent  time  has  elapsed 
so  that  the  second  argument  of  f  is  negative  at  z  >  0.  This  time  delay 
is  given  by 


6T  -  (34>/3k  )  /  (3a)/3k  ). 

z  z 


UKCLA 


When  the  field  reappears,  the  wave  packet  has  traveled  a  lateral 
distance  6L  obtained  by  setting  the  first  argument  of  f  equal  to  zero. 

The  transverse  displacement  is  given  by 

6L  -  -3(t./3k^  +  (3(i)/3k^)  6T. 

This  lateral  mode,  or  "treetop  mode,"  dominates  at  large  distances 
because  it  propagates  mainly  in  the  lossless  medium  rather  than  in  the 
jungle.  The  effect  will  be  a  large  one  only  when  the  phase  41  is  a 
rapidly  varying  function  of  either  the  frequency  or  angle.  In  our 
problem  this  occurs  at  the  critical  angle  6^  of  total  internal  re¬ 
flection.  At  0  =  0^,  6L  =  “  and  is  a  rapidly  varying  function  of  the 

difference  0  -  0  . 

c 

By  looking  in  detail  at  the  dependence  of  6L  on  0  it  is  possible 

-1/2  -3/2 

to  demonstrate  that  the  field  strength  varies  as  6r  6L  which  is 

2 

consistent  with  the  analytically  obtained  dependence  cf  1/r  . 

The  above  analysis  is  by  no  means  novel  or  original.  The  phe¬ 
nomenon  is  a  general  characteristic  of  wave  motion  and  whenever  the 
phase  of  a  wave  is  a  function  of  any  parameter  there  will  be  a  displacement 
in  the  corresponding  "canonically  conjugate"  parameter.  In  ele-  <;rical 
network  theory  the  time  delay  associated  with  a  phase  change  which 
varies  rapidly  with  frequency  is  very  familiar.  In  quantum  mechanics 
the  rapid  energy  variation  of  a  phase  shift  at  a  resonance  is  associated 
with  the  time  the  particle  spends  in  the  resonant  state.  The  lateral 
displacement  of  a  light  beam  has  been  observed  for  rays  near  the 
critical  angle,  and  this  kind  of  a  wave  is  well  known  in  seismology. 

In  our  discussion  it  was  necessary  to  make  the  assumption  that 
the  magnitude  of  the  reflection  coefficient  varied  slowly  with  angle. 

Such  an  assumption  is  correct  for  a  lossless  medium  but  is  :r'ot  valid 
if  there  are  significant  losses.  Such  losses  would  change  the  shape  of 
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the  wave  packet.  Other  effects,  such  as  a  spherical  divergence  of  the 
wave  packet,  which  change  the  shape  of  the  wave  packet  have  also  been 
ignored.  The  formal  solution  correctly  contains  all  these  effects.  As 
a  result  these  coniiideratlons  are  an  aid  in  the  understanding  of  the 
analytical  calculations  but  do  not  In  any  sense  replace  them. 

The  formal  solution  Is  originally  expressed  as  an  Integral  over 
real  a.  In  this  form  the  numerical  quadrature  described  in  Appendix  II 
Is  appropriate.  The  Integral  may  also  be  considered  as  a  contour 
Integral  In  the  complex  plane.  By  deforming  the  contour,  contributions 
from  discrete  poles  and  branch  cuts  may  be  Identified  and  separately 
considered.  This  is  discussed  in  detail  In  Appendix  III.  The  pole 
terms  correspond  to  damped  waveguide  modes.  For  sufficiently  large  r 
the  contributions  from  the  branch  cuts  correspond  to  the  lateral  waves 
as  described  above.  At  smaller  r  the  contribution  from  the  branch 
cuts  has  a  more  complicated  behavior  than  described  above  and  represents 
an  extension  of  these  simple  Ideas. 

III.  SAMPLE  CALCULATIONS  AND  PRELIMINARY  COMPARISON  WITH  EXPERIMENT 

Figure  3  Is  a  plot  of  calculated  E  (rms)  vs  range  for  a  frequency 

z 

of  0.88  Me.  For  this  low  frequency,  the  nuonerical  Integration  described 
in  Appendix  II  must  be  used  throughout  the  range  of  Interest.  Figure  3 

Illustrates  three  distinct  stages  of  the  behavior  of  E  vs  r.  For 

2  * 
r  <  10  m,  E^  decreases  exponentially.  In  this  region  the  waves  remain 

in  the  jungle  and  are  described  mathematically  by  the  pole  terms.  For 
2  3 

10  m  <  r  <  10  m,  E  decreases  as  1/r.  The  pole  terms  are  now  negligible 

compared  to  the  lateral  wave  which  has  a  1/r  dependence  at  short  ranges. 
3  4 

For  10  m  <  r  <  10  m,  the  r  dependence  of  the  lateral  wave  Is  In  a 

2 

transition  region  which  will  lead  to  the  final  dependence  of  1/r  given 

L 

by  Eq.  (40),  Appendix  III  for  r  >  10  m. 

Figure  4  Illustrates  this  transition  for  f  ••  2  Me.  The  points 
marked  with  an  x  are  the  result  of  the  Gaussian  Integration  of  the 
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Figure  4.  Sauple  Calculated  Field  for  f  «  2.0  Mc/sec 
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lateral  wave  described  in  Appendix  III.  For  this  frequency,  the  lateral 
wave  is  the  dominant  contribution  throughout  the  range  of  Interest. 

However,  it  is  best  evaluated  by  the  method  of  Appendix  II  for 

3  2 

r  <  5  X  10  m.  From  Fig.  4,  the  lateral  wave  falls  off  slower  than  1/r  , 

2  2 

then  faster  than  1/r  ,  and  finally  settles  to  a  1/r  range  dependence. 

The  remainder  of  this  section  is  concerned  with  a  preliminary 

comparison  of  the  calculated  electric  field  with  the  experimental  results 

3  4 

of  Jansky  and  Bailey.  *  The  only  case  under  consideration  at  this 
time  is  that  of  vertical  polarization  with  both  transmitter  and  receiver 
within  the  jungle.  The  frequency  range  of  this  comparison  is  from  6 
to  100  Me. 

The  distance  from  the  transmitter  is  limited  to  one  mile.  For 
ranges  larger  than  this  the  transmitting  and  receiving  antenna  were 
often  separated  by  intervening  hills,  and  the  results  do  not  directly 
apply. 


The  simple  expression  [Eq.  (40),  Appendix  III]  for  the  lateral 
wave  is  accurate  to  3  db  in  the  regime  covered  by  the  following  graphs. 
The  points  on  th'?  graphs,  however,  are  calculated  by  the  Gaussian 
integration  and  are  accurate  to  0.1%. 

For  comparison  purposes  the  calculated  electric  field  is  trans- 

4 

formed  into  basic  path  loss,  by  the  equation 

L.  -  139.0  -  20  log  E  +20  log  f 
D  rms 

where  f  is  the  frequency  in  megacycles  and  (microvolt/meter)  is 

calculated  for  a  point  dipole  radiating  one  kilowatt  into  empty  space. 

Figure  5  shows  a  series  of  curves  of  L,  vs  r  (miles)  for  various 

D 

receiver  heights  at  a  given  frequency,  f,  and  transmitter  height,  z^. 
The  straight  lines  are  calculated  values  for  the  input  parameters. 
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TRANStllSSION  FREQUENCY,  f  «  6  Mc/sec 
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SLAB  HEIGHT,  h  *  40  ft. 
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_  __  n  =  0.2  mllllmhos/m 


Figure  5,  Comparison  of  Predicted  and  Measured  Path  Loss, 

f  ■  6  Mc/sec,  z  «  21  ft 
o 


and 


h  -  40  ft 

e  «  15,  0  "10  millimho/m 
S  g 

Cj  •  1.02 


0.1  mlllimho/m  for  the  solid  line  and 


0.2  mlllimho/m  for  the  dashed  line. 


Figures  6  through  12  are  additional  graphs  for  other  combinations 
of  f  and  z^.  In  all  graphs  the  experimental  points  are  denoted  by  the 
letters  A  and  B.  These  points  denote  measurements  along  two  different 
winding  trails  (A  and  B)  at  the  Thailand  jungle  test  site.  The  scatter 
of  the  experimental  data  points  is  believed  due  to  the  lack  of  homo¬ 
geneity  of  the  jungle  and  the  variation  in  environment  of  the  antenna. 
The  environment  of  the  transmitting  antenna  can  cause  an  azimuthal 
asymmetry  in  the  radiation  patterns.  Also,  the  jungle  character  of 
the  ray  path  at  each  azimuth  may  be  different.  Thus,  there  will  be  an 
azimuthal  asymmetry  in  the  electric  field.  Since  the  trails  are  not 
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Figure  8.  Comparison  of  Predicted  and  Measured  Path  Loss, 
f  ■  25.5  Mc/sec,  =  40  ft 
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Figure  10.  Comparison  of  Predicted  and  Measured  Path  Loss 
f  “  50  Mc/sec,  z  *  40  ft 


Since  the  calculated  electric  field  depends  most  critically  on 
the  jungle  conductivity,  this  was  the  only  quantity  varied  in  the 
calculations.  The  curves  for  the  two  values  of  conductivity  generally 
bound  the  experimental  points  for  all  values  of  receiver  and  trans¬ 
mitter  height  and  frequency  with  the  exception  of  the  f  ■  25.5  Me, 

l^ 

■  21  ft  case  (Fig.  7).  Jansky  and  Bailey,  on  pp.  4.6  -  4.8  of 
Report  6,  call  attention  to  an  apparent  inconsistency  in  their  data 
at  25.5  Me  and  low  transmitter  height  and  observe  that  a  10  db  decrease 
in  their  path  loss  would  make  their  results  more  reasonable.  This 
10  db  decrease  in  measured  path  loss  will  bring  the  experimental 
points  into  the  neighborhood  of  the  calculated  values. 
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DIFFERENCE  BETWEEN  PREDICTED  AND  fiEASUREO  BASIC  PATH  LOSS  VALUES  (db) 


The  difference  between  the  calculated  and  measured  path  loss  is 
plotted  statistically  in  Fig.  13  for  the  two  conductivities.  The 
points  of  Fig.  7  were  not  included.  These  plots  show  the  lower 
conductivity  to  predict  a  path  loss  which  is  too  low  by  4  db  and  the 
higher  conductivity  to  predict  a  loss  high  by  5  db.  The  important 
result  is  the  standard  deviation  of  6  db.  This  Illustrates  the  utility 
of  the  slab  model.  The  calculation  has  a  68Z  chance  of  being  within 
6  db  of  the  measurement  for  all  antenna  heights,  ranges  and  frequencies 
between  6  and  100  Me  when  the  Jungle  conductivity  is  taken  as  approxi¬ 
mately  0.15  mllllmho/m. 


PERCENT  OF  SAMPLE  EXCEEDING  ORDINATE  VALUE 


Figure  13.  Distributions  of  Differences  Between  Predicted 
and  Measured  Path  Loss 


IV.  CONCLUSIONS 


A  simple  slab  model  of  the  Jungle  is  capable  of  explaining  the 
experimental  propagation  data  over  a  wide  range  of  parameters.  The 
calculations  performed  so  far  have  considered  only  vertical  polarization 
and  both  antennas  in  the  jungle.  Me  are  continuing  the  study  to  In¬ 
clude  horizontal  polarization  and  cases  where  either  or  both  of  the 
antennas  are  above  the  jungle.  A  comparison  with  the  experimental 
results  of  Jansky  and  Bailey  and  other  experimental  groups,  where 
available,  will  be  made  and  the  regimes  of  validity  of  the  slab  model 
will  be  determined. 
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APPENDIX  I 

DERIVATION  OF  THE  FIELD  EXPRESSIONS 

The  integral  representation  of  the  fields  can  be  obtained  simply. 

For  the  sake  of  completeness  we  include  an  outline  of  the  derivation 
following  that  given  by  Walt. 

The  free-field  equations  are 

V  X  E  -  i  -  H  (8) 

c 

0X1/  i(»ig  .  Atto  .  _  iw  2„ 

V  X  H  »  (-  -7-  +  )  E  «  -  n  E 

c  c  c 

where  u  ■  2irf  •  cK^  is  the  angular  frequency  of  the  radiation.  The  fields 
can  be  expressed  in  terms  of  the  Hertz  vector  n 

E  »  V  X  (V  X  n)  (9) 

H  -  -  iK^n^V  X  n  (10) 

if  n  satisfies  the  equation 

v^n  +  n^K  -  0.  (11) 

O 

Because  of  the  cylindrical  symmetry  of  the  problem  (with  vertical  polari¬ 
zation),  n  is  in  the  z-directlon. 

The  solutions  of  Eq.  (11)  which  are  regular  at  the  origin  and 
Independent  of  azimuth  angle  are 
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n  -  J^  (aR)  e“  (12) 

2  2  2 

where  R  Is  measured  In  units  of  c/cu  and  x  ^  a  -  n  ,  and  the  real 
part  of  X  Is  greater  than  zero. 


The  solution  to  our  problem  satisfies  the  free-fleld  equations 
everywhere  except  at  the  antenna  where  It  has  a  singularity  character¬ 
istic  of  a  point  dipole.  By  superposition  the  solution  has  the  form 


»  X  (H-Z) 

^  s_(a)  J  (aR)  e  ada 

•'o 


a  o 


(13) 


vXE 


-X  X  X  (H-Z)  -X  (H-Z)]  oJ  (aR) 

a  e  ^  +  a,  (a)  e  ^  +  b .  (a)  e  ^  I -  da 

o  J  3  J 


j 


(14) 


r  -X  (H-Z) 

n  -  I  b  (a)  e  ®  J  (aR)a  da 
g  I  g  o 

*'o 


(15) 


where  the  coefficient  a^  Is  a  measure  of  the  dipole  moments  of  the  trans¬ 
mitting  antenna  and  the  other  symbols  have  the  same  significance  as  In 
Section  II. 

By  using  Eqs.  (9)  and  (10)  the  field  components  can  be  determined 

as  linear  functions  of  the  four  unknown  functions  a  (a) ,  a . (a) ,  b , (a) ,  b  (a) . 

^  3  3  g 

The  requirement  that  the  tangential  components  of  the  electrical  and 
magnetic  field  be  continuous  at  the  boundary  leads  to  four  algebraic 
equations  vdilch  determine  these  coefficients.  The  algebra  Is  straight¬ 
forward  and  leads  to 

00 

E  -  2P  I  a^  J  (oR)F(aR)da  (16) 

2  10 

•'o 
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where  P  and  F(aR)  are  as  given  in  Eq.  (1)  in  the  text.  It  is  convenient 
to  have  the  range  of  a  extend  from  -«>  to  +».  To  do  this  we  note  that 
F(-a)  *  -F(o)  and 

J^(aR)  -  j  (aR)  -  H^^^^(-aR)]. 

With  these  substitutions  we  obtain  Eq.  (1)  of  the  text. 


APPENDIX  II 

EXPLICIT  EVALUATION  OF  THE  FIELD  INTEGRAL 

The  explicit  evaluation  of  the  field  Integral  given  by  Eq.  (1)  of 
Section  II  cannot  be  achieved  by  conventional  quadrature  techniques 
because  of  the  rapid  oscillation  of  the  Bessel  function  for  large  aR. 

Thus  even  though  the  function  F(o)  converges  in  general  within  modest 
ranges  of  a,  most  numerical  Integration  procedures  would  lose  all  signif¬ 
icance  whenever  R  became  appreciable  because  of  the  cancellation  of  the 
individual  quadrature  terms  with  one  another.  For  any  well-behaved  F(a) 
the  Integral  Is  easily  shown,  to  tend  to  zero  as  R  tends  to  infinity.  Via 
numerical  quadrature,  this  would  be  manifest  by  a  large  number  of  terms 
of  essentially  equal  magnitude  (though  different  sign)  being  added  together 
to  yield  a  residual  which  may  be  made  as  S!.>all  as  desired  by  letting  R 
become  large.  Irrespective  of  the  number  of  significant  figures  a  given 
digital  computer  is  capable  of  handling,  there  will  always  exist  an  R 
beyond  which  a  mechanical  quadrature  will  result  in  a  total  loss  of 
significance.  There  is,  however,  a  means  of  numerically  evaluating  the 
integral  in  question  without  meaningful  loss  of  accuracy  for  any  value  of 
R  so  long  as  computing  time  is  no  obstacle.  In  other  words,  a  numerical 
attack  upon  the  integral  is  always  possible,  but  becomes  more  time- 
consuming  as  R  increases.  The  broad  range  of  validity  for  the  analytic 
methods  of  Section  III,  however,  insures  that  these  contingencies  do  not 
arise.  Of  particular  note  concerning  the  numerical  evaluation  described 
herein  is  its  exact  agreement  with  the  results  of  Section  III  in  the 
regime  beyond  which  the  latter  approach  was  shown  to  be  accurate.  Thus 
two  essentially  independent  evaluations  of  the  integrals  concerned  yield 
identical  results,  thereby  providing  a  self-consistency  check  rarely 
available  in  similar  work.  In  regimes  where  the  methods  of  Section  III 
are  no  longer  accurate,  those  described  in  this  appendix  ^re  used;  and 
conversely,  where  these  numerical  quadratures  become  too  time-consuming, 
the  analytic  methods  of  Section  III  are  of  sufficient  accuracy  to  be  used. 
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A  mechanical  quadrature  of  the  field  Integral  Is  achieved  by  first 
noting  that  the  source  term  contribution  may  be  evaluated  explicitly  • 
Thus 

00 

E,  -  E  +  A  I  If.  (a)  +  f-(a)]  J  (aR)o^do  (17) 

Z  Z  f  1  l  o 

'  •) 


where 


(18) 


and 


Q  -  Vr^  +  (Z-Z  )^ 


(19) 


with  f^(a)  and  f2(a)  given  In  appendix  I. 


2  2 
When  a  becomes  somewhat  larger  than  the  irodulus  of  n^  ,  the  dominant 


behavior  of  f^(a)  Is  controlled  by  the  exponential  factor 


-■‘j’'  -aY 
e  e 


(20) 


whereas  £2(0.)  Is  dominated  by  the  factor 
-X, (2H-Y) 


j 


-a(2H-Y) 


(21) 


Thus  given  h,  z,  z^,  and  f  (See  Sec.  II  for  nomenclature)  an  upper  limit 

of  a  (called  a  )  may  be  chosen  beyond  which  the  Integrand  of  Eq.  (17) 
mdx 

contributes  a  negligible  amount.  This  a  will  be  sufficient  for  a 

DldX 
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complete  range  of  R.  If  R  is  fixed,  on  the  other  hand  (i.e.,  the  genera¬ 
tion  of  height  gain  curves  is  required) ,  must  be  varied  in  such  a  way 

that  Ojugjj  Y  is  greater  than  a  fixed  constant.  Since  Y  and  H  are  directly 

proportional  to  the  propagation  frequency,  a  is  easily  frequency  scaled. 

indx 

All  that  remains,  therefore,  is  the  evaluation  of 
g 

/max 

[fj^(g)  +  f2(o)]  J^(gR)g^dg  (22) 

o 


where  fj^(g)  and  f2(oi)  are  complex  functions.  Separating  each  into  real 
and  imaginary  parts,  the  secondary  field  may  still  be  written  in  terms  of 
finite  Integrals  along  the  real  g-axis  as 


sec 


r  o. 

/max 

8j^(a) 

L  o 


J  (gR)g' 
o 


dg  +  i 


U 

/ 


max 


g2(“) 


J  (gR)g' 
o 


da 


(23) 


where  gj^(cx)  -  Retfj^(g)  +  f2(a)]  and  *  Im  [fj^(g)  +  f2(a)]. 

The  severe  oscillations  of  J  (gR)  do  not  occur  until  gR  is  well 

o 

beyond  unity.  Indeed  at  this  point  J^(gR)  can  be  accurately  replaced  by 
the  first  two  terms  of  its  asymptotic  expansion,  viz. 

vis  -.."tf)].  (24, 


This  expansion  will  provide  an  approximation  to  J  (oR)  better  than  one  part 

L  ° 

in  10  when  gR  >  30.  Choosing  6  =  30/R  the  integrals  of  Eq.  (23)  may  be 
divided  into  two  parts,  viz. 


6  5max 

/  -/ 

•'o  6 


dg  , 


(25) 
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with  the  asymptotic  representation  of  J^(aR)  being  used  In  the  last  inte¬ 
gral.  The  Integrals  from  0  to  3  may  be  very  accurately  generated  using 
standard  Gaussian  quadratures  whereas  the  Integrals  from  3  to  are 
most  amenable  to  solution  by  Filon's  method.^  This  latter  method,  des¬ 
cribed  in  considerable  detail  in  references  9  and  10,  was  developed  to 
handle  mechanical  quadratures  of  the  type 


f(x)  cos  Rx  dx 


Upon  substituting  (24)  into  (23)  for  a  >  3,  the  resulting  integrals 
may  be  cast  into  one  of  the  forms 


/max 


(x)  cos  (Rx  -ir/4)dx  or 


F2(x)  sln(Rx  -  iT/4)dx  (27) 


The  Interval  from  3  to  a  is  then  divided  into  2n  equal  intervals  of 

max  ^ 

size  6  *  6)/2n  and  the  values  of  the  indepen'^ent  variable  o  at 

these  points  are  indicated  by  x^,  Xj^,  ...,  X2^.  Introducing  the  variable 
e  ■  dR,  one  obtains 


F(a)  cos(aR  -  TT/4)da 


‘  l”2n  *‘"‘“*211  ■  ■  “o  -  f)  t 


+  q(e)C2„  +  r(e)C2^.J  +  r„ 
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F(o)  sln(aR  -  Tr/4)da 


«[p(e)  j",  cos<Rx^  -  i)  -  cosdiXj^  -  J)( 


+  1(6)  S2„  +  F(e)S2„.J+ 


where  the  generalized  Slmpaon  coefficients  p,  q,  and  r  are  given  by 


P(e)  .  i  4 

26  6 


q(6)  -  2 


cos  6  sin  29 


! 

I  ' 


1(0)  -  4 


'sin  9  cos  0 


)  • 


In  addition,  one  has 


■  F2^  cos  (Rx2^  -  Tr/4) 


~  2  P2n  ^*^2n  "  ^o  ‘'°®  ^^o  ' 


‘^2n-l  ■  2  ''21-1  ‘“*21-1  ■ 


hr,  ■  2  “21  ‘“*21  ■  "  2  ‘“20  ‘“*2n  ' 


+  F  sin  (Rx  -  7r/4)l 
o  o 
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iiAioi  lec^inrii 


s 


2n-l 


n 

X!)  ^21-1  ^^2i-l  “ 

i-1 


(36) 


As  la  evident  from  expressions  (37)  and  (38) ,  the  error  Is  propor¬ 
tional  to  R  times  the  Interval  size  to  the  fifth  power.  Thus  for  very 
large  R  the  Interval  size  h  must  be  decreased  appropriately  to  maintain 
accuracy.  This  In  turn  results  In  Increasingly  long  periods  of  computa¬ 
tion.  The  true  practical  usefulness  of  Fllon's  method  Is  therefore  In 
some  doubt.  Although  postulated  to  be  a  panacea  for  large  R,  the  method 
Is  actually  most  useful  from  the  practical  point  of  view  for  small  6,  l.e., 
for  hR  small.  With  a  computer,  however,  the  method  does  yield  a  convergent 
solution  as  the  interval  size  h  tends  to  zero.  This  should  be  contrasted 
to  any  other  mechanical  quadrature  scheme  which,  irrespective  of  the 
smallness  of  the  Integration  step  size,  will  in  general  not  converge  for 
any  reasonably  large  valoe  of  R. 


APPENDIX  III 
ANALYTIC  SOLUTIONS 

In  order  to  obtain  an  asymptotic  expansion  of  the  fields  we  first 
note  that  has  the  asymptotic  expansion 

“o”  <“«  -  vis  (i  ks  “  (^)) 

valid  for  |aR|  >>  1.  This  may  then  be  inserted  into  Eq.  (1)  of  the  text 
and  the  resulting  Integral  can  be  evaluated  as  a  contour  integral. 

The  Integrand  vanishes  on  an  infinite  semi-circle  In  the  upper 
half  a-plane  because  of  the  term  in  .  The  original  contour  which 
lay  along  the  real  axis  is  therefore  deformed  to  this  semi-circle.  One 
is  left  with  Integrals  along  two  branch  cuts  and  contributions  from  poles 
that  were  crossed  (Fig.  14). 


The  poles  correspond  to  transmission  modes  through  the  Jungle  which 

are  exponentially  damped.  The  cut  at  x  »  0  (a=  n  )  corresponds  to  the 

E  E 

lateral  wave  at  the  jungle-ground  boundary.  This  wave  is  highly  attenu¬ 
ated  because  of  the  high  conductivity  of  the  ground  and  is  negligible. 

The  cut  at  X  »  0  (a  >  1)  corresponds  to  the  lateral  wave  at  the  jungle- 
air  boundary.  This  wave  is  not  exponentially  attenuated  and  is  therefore 
the  dominant  contribution  after  one  damping  length  of  the  pole  terms. 


There  is  no  branch  point  at  x. 


n^ )  since  F  is  even  in  x^ . 


The  poles  are  the  normal  modes  of  propagation  through  a  waveguide 
and  would  be  the  dominant  contribution  if  the  waveguide  were  not  a  lossy 
jungle.  In  this  case  these  pole  contributions  are  negligible  for 
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Figure  14.  Contours  in  the  Complex  o-Plane 
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distances  greater  than  a  jungle  wavelength.  The  cut  at  a  >  1  corresponds 
to  that  part  of  the  sovrce  radiation  which  is  transmitted  from  the  jungle 

o 

to  the  air  above  and  then  propagates  through  the  air.  This  is  analogous 
to  the  usual  ground  or  surface  wave  of  radio  propagation  except  in  this 
case  the  lower  boundary  is  jungle  rather  than  ground. 


The  contribution  from  the  cut  at  a  ■  1  is  evaluated  by  transforming 

2 

to  the  variable  s  by  the  transformation  a«l  +  l8  ;-<»£s£9’.  This 
transformation  describes  the  contour  of  the  a-plane  which  runs  from  a  «  1 
+  i<»  to  a  *  1  on  the  left  of  the  cut  as  the  negative  real  s  axis  and  the 
contour  in  the  a-plane  which  runs  from  a  -  1  to  o  »  1  +  i*  on  the  right 
side  of  the  cut  as  the  positive  real  s  axis. 

By  taking 

X  ■  V -  1  ■  sd  +  1)  V  1  +  is^/Z 
a 


the  correct  sign  for  x  on  each  side  of  the  cut  obtains. 
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With  this  transformation  the  branch  cut  contribution  is 

on 

2 


JT  iTr/A  iR 
® 


I  ds(l  +  is 
•^00 


2v3/2  "Rs  f  \\/i  • 

)  8  e  F(o(s))(l  + 


81R(l+is  ) 


+  . . .) 


(39) 


2  ro 

The  functions  F(s)  and  (1  +  is  )  where  m  =  -1  and  5/2  are  expanded 
in  Taylor  series  about  s  »  0.  Only  terms  of  the  Integrand  that  are  even 
in  s  survive.  The  result  is 


where  the  subscript  "o"  means  the  derivatives  are  evaluated  at  s  <■  0.  The 
expression  for  is  an  asymptotic  expansion  in  R.  The  ^  terms  in  the 
bracket  provide  an  estimate  of  the  error  Involved  when  using  the  first 
term  alone.  The  result  is  then 

,PF  ’  , 

E  ■ 
rms  p2 


9aO^^Vpower 

VTvf  ln^-l|r^ 


K  (Z,  Z^,  H)  pvolt/m 


(AO) 


where 
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This  expression  is  too  complicated  to  be  useful  If  a  quick  estimate 
of  the  error  Involved  In  using  Eq.  (40)  Is  desired.  Comparison  with  a 
numerical  evaluation  of  the  branch-cut  Integral,  Eq.  (39),  shows  that 
the  discrepancy  between  Eq.  (40)  and  the  exact  evaluation  agrees  with 
the  estimated  error  only  at  the  high  frequencies  (f  ^  25.5  Me).  At 
f  <  25.5  Me,  the  discrepancy  Is  greater  than  the  estimated  error  by  a 
factor  of  three  In  some  cases.  It  would  therefore  be  desirable  to 
replace  the  error  expression,  Eq.  (41),  by  a  simpler  and  more  accurate 
error  bound.  Numerical  work  shows  that  the  expression  for  the  error  may 
be  approximated  by  the  following  formula: 


* 

3  (H  -  Y/2) 

R(nj^  -  1) 

+ 

(42) 


for  frequencies  in  the  range  6  <  f  <  100  Me.  However,  this  expression  Is 
overly  pessimistic  at  the  high  frequencies  and  Is  too  small  for  f  <  6  Me. 
An  accurate  error  bound  at  all  frequencies  can  be  obtained  by  simply 
using  three  times  E  (expression  41).  More  extensive  numerical  comparisons 
will  be  necescary  In  order  to  determine  simpler  expressions  for  accurate 
error  bounds  In  all  frequency  regions. 

Both  expressions  (41)  and  (42)  show  the  error  to  be  Inversely 
proportional  to  r.  The  error  decreases  as  the  range  Increases.  Thus  the 
long-range  behavior  of  the  electric  field  Is  accurately  given  by 
expression  (40) . 

When  E  =  1,  expression  (40)  Is  no  longer  accurate  but  the  integrand 
of  Eq.  (39)  is  still  well-behaved  so  that  a  Gaussian  Integration  of 
Eq.  (39)  is  easily  effected.  In  this  regime,  the  branch  cut  is  still 
the  dominant  contribution.  However,  it  is  evaluated  numerically  to  an 
accuracy  of  0.1%  rather  than  by  the  asymptotic  expansion. 
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For  E  >  1,  the  Gaussian  integration  falls  to  converge.  In  this 

regime  a  numerical  integration  of  the  original  integral  is  used.  The 

numerical  integration  is  accomplished  by  a  machine  code,  described  in 

Appendix  II,  which  evaluates  the  total  field,  that  is,  the  contribution 

from  the  cuts  and  all  poles.  We  thus  are  in  a  position  to  make  an 

accurate  numerical  evaluation  of  E  for  all  values  of  the  input  parameters. 
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A  theoretical  determination  of  the  path  loss  in  radio  propagation  through 
jungle  Is  obtained  by  considering  the  jungle  as  a  homogeneous  conducting 
dielectric  slab  on  a  flat  earth.  An  exact  Integral  is  obtained  for  the 
vertical  component  of  the  electric  field  within  the  jungle  due  to  a  vertical 
electric  dipole  within  the  jungle.  An  analytic  evaluation  of  the  Integral 
leads  to  an  approximate  expression  for  the  field  and  an  estimate  of  the 
error.  When  this  accuracy  is  insufficient,  a  numerical  evaluation  of  the 
Integral  is  performed.  The  combination  of  analytic  and  numerical  techniques 
leads  to  the  evaluation  of  the  electric  field  and  thereby  path  loss  to  any 
accuracy  desired.  Experimental  measurements  of  path  loss  in  a  Thailand 
jungle  are  compared  with  the  calculations.  The  results  agree  within  a 
standard  deviation  of  6  db  when  the  jungle  conductivity  is  taken  as  0.15 
milllmho/m. 
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